Double neutron/proton ratio of nucleon emissions in isotopic reaction systems as a 

robust probe of nuclear symmetry energy 
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The double neutron/proton ratio of nucleon emissions taken from two reaction systems using 
four isotopes of the same element, namely, the neutron/proton ratio in the neutron-rich system 
over that in the more symmetric system, has the advantage of reducing systematically the influence 
of the Coulomb force and the normally poor efficiencies of detecting low energy neutrons. The 
double ratio thus suffers less systematic errors. Within the IBUU04 transport model the double 
neutron/proton ratio is shown to have about the same sensitivity to the density dependence of 
nuclear symmetry energy as the single neutron/proton ratio in the neutron-rich system involved. 
The double neutron/proton ratio is therefore more useful for further constraining the symmetry 
energy of neutron-rich matter. 

PACS numbers: 25.70.-z, 24.10.Lx 



I. INTRODUCTION 



The density dependence of nuclear symmetry energy 
Esym{p) is still poorly known but very important for both 
nuclear physics and astrophysics P, 12, 0,0, IE IE| • Heavy- 
ion reactions induced by neutron-rich nuclei provide a 
unique opportunity to constrain the symmetry energy in 
a broad density range. A number of potentially useful 
probes of the symmetry energy have been proposed in the 
literature mostly based on transport model simulations. 
Experimental data from several dedicated experiments 
just started emerging. Some comparisons between the 
available experimental data and transport model calcu- 
lations have been carried out recently. These studies have 
allowed us to place important constraints on the density 
dependence of symmetry energy. For instance, by using 
the free-space experimental nucleon-nucleon (NN) cross 
sections within the transport model IBUU04 [3, a sym- 
metry energy of -Esym(/3) ~ 31.6(p/ poY'^ for densities less 
than 1.2pQwas extracted from the MSU data on isospin 
diffusion 0, |||. While using in-medium NN cross sec- 
tions calculated within an effective-mass scaling approach 
[13, a symmetry energy of Esymip) « 31.6(p/po)" ''^ 
was found most acceptable in comparison with both the 
MSU isospin diffusion data and the presently acceptable 
neutron-skin thickness in ^'^^Pb [Tol [lll|. The currently 
existing isospin diffusion data alone can not distinguish 
the above two forms of the symmetry energy within the 
experimental error bars Thus complementary ob- 

servables sensitive to the £'sym(p), more desirably, studies 
on correlations of several such observables, are still very 
much needed to further constrain the symmetry energy. 

Because the symmetry potentials have opposite signs 
for neutrons and protons and the fact that the symmetry 
potentials are generally smaller compared to the isoscalar 
potential at the same density, most of the observables 



proposed so far use differences or ratios of isospin multi- 
plets of baryons, mirror nuclei and mesons, such as, the 
neutron/proton ratio of nucleon emissions J^], neutron- 
proton differential flow |l3|. neutron-proton correlation 
function 111, t/^He '^'-I^^ El 111 111 113, 

I]-/E+ liirand K°/K+ ratios H^, etc. Among these 
observables, the neutron/proton ratio of nucleon emis- 
sions has probably the highest sensitivity to the sym- 
metry energy. This is because symmetry potentials act 
directly on nucleons and normally nucleon emissions are 
rather abundant in typical heavy-ion reactions. However, 
it is very challenging to measure some of these observ- 
ables, especially those involving neutrons. The measure- 
ment of neutrons, especially the low energy ones, always 
suffers from low detection efficiencies even for the most 
advanced neutron detectors. Therefore, observables in- 
volving neutrons normally have large systematic errors. 
Moreover, for essentially all of these observables, the 
Coulomb force on charged particles plays an important 
role, ft sometimes competes strongly with the symme- 
try potentials. One has to disentangle carefully effects of 
the symmetry potentials from those due to the Coulomb 
potentials. It is thus very desirable to find experimental 
observables which can reduce the influence of both the 
Coulomb force and the systematic errors associated with 
neutrons. The double neutron/proton ratio of nucleon 
emissions taken from two reaction systems using four iso- 
topes of the same element, namely, the neutron/proton 
ratio in the neutron-rich system over that in the more 
symmetric system, was recently proposed by Lynch et 
al. 2s| as a candidate of such an observable. They have 
actually measured the double neutron/proton ratio in 
central reactions of ^^"'Sn-l-^^^'Sn and ^^^Sn-f^^^Sn at a 
beam energy of 50 MeV/nucleon at the National Super- 
conducting Cyclotron Laboratory '2^. While the experi- 
mental data are currently being finalized, we report here 
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transport model analyses of the double neutron/proton 
rations. Besides the above two reactions we also study 
the double neutron/proton ratio in ^^^Sn+^^'^Sn and 
^^^Sn+^^^Sn reactions at 400 MeV/nucleon. It is shown 
that the double neutron/proton ratio has about the same 
sensitivity to the density dependence of symmetry en- 
ergy as the corresponding single ratio in the respective 
neutron-rich system involved. Given the advantages of 
measuring the double neutron/proton ratios over the sin- 
gle ones, the study of double neutron/proton ratios will 
be more useful for further constraining the symmetry en- 
ergy of neutron-rich matter. 



II. A SUMMARY OF THE IBUU04 
TRANSPORT MODEL 

Our study is carried out using the IBUU04 version of 
an isospin and momentum dependent transport model 
for nuclear reactions induced by neutron-rich nuclei 0. 
For completeness and consistency we outline here a few 
major features most relevant to the present study. More 
details of the model can be found in Refs. The 
single nucleon potential is one of the most important in- 
puts to all transport models for nuclear reactions. In the 
IBUU04 transport model, we use a single nucleon po- 
tential derived within the Hartree-Fock approach using a 
modified Gogny effective interaction (MDI) i.e.. 
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Here 5 = (pn — Pp) /p is the isospin asymmetry of the nu- 
clear medium. In the above r = 1/2 (—1/2) for neutrons 
(protons) and t ^ t' ] a = fr{r,p) is the phase space 
distribution function at coordinate r and momentum p. 
The parameters Au{x)^ Ai[x), B ,Ct-.t,Ct-.t' and A were 
obtained by fitting the momentum-dependence of the 
U{p,S,p,T,x) to that predicted by the Gogny Hartree- 
Fock and/or the Brueckner-Hartree-Fock (BHF) calcula- 
tions 25] , the saturation properties of symmetric nuclear 
matter and the symmetry energy of about 30 MeV at 
normal nuclear matter density po = 0.16 fm~'^ |2^. The 
incompressibility Kq of symmetric nuclear matter at po 
is set to be 211 MeV consistent with the latest conclu- 
sion from studying giant resonances [2^ [23, The 
parameters Au{x) and Ai{x) depend on the x parameter 
according to 



Au{x) = -95.98- 
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Ai{x) = -120.57- 
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The parameter x can be adjusted to mimic predictions on 
the density dependence of symmetry energy Egymip) by 
microscopic and/or phenomenological many-body theo- 
ries. Shown in Fig. is the density dependence of the 
symmetry energy for a; = and —1. The recent anal- 
yses of the MSU isospin diffusion data have allowed us 
to constrain the x parameter to be between these two 
values for densities less than about 1.2po 10]- The cor- 
responding symmetry energy can be parameterized as 
Esyn.{p) « 31.6(p/po)'-' and E,y^{p) « 31.6(p/po)°-^' 
for X = — 1 and a; = 0, respectively. The main pur- 
pose of this work is to investigate whether the double 
neutron/proton ratio can help further narrow down the 
uncertainty of the symmetry energy. 




FIG. 1: (Color online) Symmetry energy as a function of den- 
sity for the MDI interaction with x = and —1. 

The last two terms in Eq. contain the momentum- 
dependence of the single-particle potential. The momen- 
tum dependence of the symmetry potential stems from 
the different interaction strength parameters Cr^r' and 
Cr,T for a nucleon of isospin t interacting, respectively, 
with unlike and like nucleons in the background fields. 
More specifically, we use CunUke = —103.4 MeV and 
Ciike = -11.7 MeV. With these parameters, the isoscalar 



potential estimated from (C/„ 



Uproton)/'^ agrees 



reasonably well with predictions from the variational 
many-body theory the BHF approach plIMIp 
including three-body forces and the Dirac-Brueckner- 
Hartree-Fock (DBHF) calculations [s^l in broad ranges of 
density and momentum. For the MDI potential we used 
here, the neutron-proton effective mass splitting due to 
the momentum dependence of the symmetry potential is 
positive 01 • This is consistent with predictions of both 
the BHF and DBHF models 25, 30, 3L ,32,, 33J. 

The IBUU04 model can use either the free-space ex- 
perimental NN cross sections [s^l or the in-medium NN 
cross sections calculated using an effective-mass scaling 
model consistent with the single particle potential used 
|lOj |. In the present work the in- medium NN cross sec- 
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tions are used. 
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FIG. 2: (Color online) Time evolution of the neutron/proton 
ratio of free nucleons as a function of kinetic energy for the 
reaction of ^^^Sn+^^^Sn (lower windows) and "^Sn+^^^Sn 
(upper windows) at 50 MeV/nucleon and an impact parame- 
ter of 2 fm, respectively. The calculations are done with the 
MDI interaction of a; = (filled square) and x — —1 (open 
square) . 



III. RESULTS AND DISCUSSIONS 

We present and discuss here the double neu- 
tron/proton ratios of free nucleons in central and mid- 
central reactions using three Sn isotopes at a beam energy 
of 50 MeV/nucleon and 400 MeV/nucleon, respectively. 
Free nucleons are identified as those having local baryon 
densities less than po/8. For a comparison and as an ex- 
ample, we first study the single neutron/proton ratio of 
free nucleons for ^^''Sn-|--^^''Sn and ^^^Sn-|-^^^Sn reactions 
at a beam energy of 50 MeV/nucleon and an impact pa- 
rameter of 2 fm, respectively. Shown in Fig. |21 are the 
time evolutions of the single neutron/proton ratios ver- 
sus the nucleon kinetic energy in the c.m.s frame of the 
respective reaction. It is seen that the neutron/proton ra- 
tio becomes stable after about 100 fm/c. As one expects, 
the neutron/proton ratio in the neutron-richer system is 
more sensitive to the symmetry energy, especially for fast 
nucleons. With the softer symmetry energy of a: = 0, the 
symmetry energy and the magnitude of the symmetry 
potential are higher at sub-saturation densities compared 
to the case with x — —1. At supra-saturation densities, 
however, it is just the opposite as shown in Fig. ^ In 
the above two reactions the maximum density reached 
is about 1.2/30 9J. One thus expects to see a higher neu- 
tron/proton ratio of free nucleons with the softer sym- 
metry energy of a; = due to the stronger repulsive 
(attractive) symmetry potential for neutrons (protons). 
For the more symmetric system ^^^Sn-f ^-"^^Sn, effects of 
the symmetry energy are negligible because of the small 



isospin asymmetry in the system. The rise of the neu- 
tron/proton ratio at low energies in both systems is due 
to the Coulomb force which pushes protons away from 
the center of mass of the reaction. These features are 
all consistent with those found in an earlier study using 
a momentum- independent transport model p^ . Unlike 
the results in Ref. T^l, however, the observed symmetry 
energy effect is only about 10% to 15% even for the most 
energetic nucleons in the ^^'*Sn-|-^^''Sn reaction. This is 
understandable since the symmetry energy used here has 
already been severely constrained by the recent isospin 
diffusion data. In ref.|l2|. however, a much wider un- 
certainty range between approximately 30(p//9o)*''^ and 
30(p/po)^ was used for the symmetry energy. 
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FIG. 3: (Color online) The double neutron/proton ratio of 
free nucleons taken from the reactions of ^■^*Sn-|-^'^*Sn and 
^^^Sn+^^'^Sn at 50 MeV/nucleon and an impact parameter of 
2 fm (upper window) and 6 fm (lower window), respectively. 

We now turn to the double neutron/proton ratio of free 
nucleons. Shown in Fig.|21are the double neutron/proton 
ratios calculated for the i24Sn-hi24Sn and "2311-^1123^^ 
reactions at a beam energy of 50 MeV/nucleon and an im- 
pact parameter of 2 fm (upper window) and 6 fm (lower 
window), respectively. As a reference, a straight line at 
74/62 corresponding to the double neutron/proton ratio 
of the entrance channel is also drawn. Below the pion 
production threshold, statistically, one would expect the 
double neutron/proton ratio of nucleon emissions to be 
a constant close to this value neglecting effects due to 
both the Coulomb and symmetry potentials. Indeed, the 
observed double neutron/proton ratios especially at low 
kinetic energies with a; = — 1 at both impact parame- 
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ters are almost constants just slightly above the straight 
line at 74/62. The fact that the double neutron/proton 
ratio is slightly higher than 74/62 even with x = — 1 
(which corresponds to a weaker symmetry potential at 
sub-saturation densities compared to the case with x=0) 
is due to the appreciable repulsive/ attractive symmetry 
potential on neutrons/protons in the ^^^Sn+^^^Sn reac- 
tion. For the double ratios in the two reactions involv- 
ing isotopes of the same element one expects that the 
Coulomb effects are largely cancelled out. More ener- 
getic nucleons have gone through denser regions of the 
reactions, effects of the symmetry potentials on them are 
thus higher especially in the case with x = 0. There- 
fore, the double neutron/proton ratios increase when the 
X parameter is changed from = — Itox^O, especially 
for energetic nucleons. 

At both impact parameters, effects of the symmetry 
energy are about 10% — 15% changing from a; = — 1 
to cc = for the energetic nucleons which are mostly 
from pre-equilibrium emissions. The observed sensitivity 
to the symmetry energy is about the same as the single 
neutron/proton ratio shown in Fig. 2. It should also be 
mentioned that since the neutron/proton ratio at kinetic 
energies less than about 50 MeV is rather insensitive to 
the symmetry energy in the reactions at a beam energy 
of 50 MeV/nucleon, neutron detectors with a threshold 
energy of 50 MeV is sufficient for the study discussed 
here. However, as we shall discuss in the following, with 
reactions at beam energies above the pion production 
threshold even the low energy neutrons are useful. 
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FIG. 4: (Color online) The double neutron/proton ratio of 
free nucleons taken from the reactions of ^^^Sn+^^'^Sn and 
^^^Sn+^^^Sn at 400 MeV/nucleon and an impact parameter 
of 1 fm (left window) and 5 fm (right window), respectively. 

Moving to beam energies above the pion production 
threshold, the reference line at 74/62 is no longer use- 
ful. In fact, the 7r~/7r+ ratio itself is a promising probe 
of the symmetry energy at high densities flTl Il3. l20j . 
Although the symmetry energy at supra-saturation den- 
sities is currently not constrained by any experimental 



data, for this study, we keep using the two parameters 
a; = and a; = —1 to be consistent with the symmetry 
energy used for sub-saturation densities. Shown in Fig. 
0] are the double neutron/proton ratios from the reac- 
tions of ^•^^Sn+^^'*Sn and ^^^Sn+^^^Sn at a beam energy 
of 400 MeV/nucleon and an impact parameter of 1 fm 
(left window) and 5 fm (right window), respectively. At 
both impact parameters, effects of the symmetry energy 
are about 5% — 10% changing from the case with x = 
to X — —1. One notices here that the low energy nu- 
cleons are having the largest sensitivity to the variation 
of the symmetry energy for such high energy heavy-ion 
collisions. In fact, the neutron/proton ratio of midrapid- 
ity nucleons which have gone through the high density 
phase of the reaction are known to be most sensitive to 
the symmetry energy |35j . Compared to the results at 
the beam energy of 50 MeV/nucleon, it is interesting to 
see a clear turnover in the dependence of the double neu- 
tron/proton ratio on the x parameter, namely the dou- 
ble ratio is lower at 50 MeV/nucleon but higher at 400 
MeV/nucleon with a; = — 1 than that with a; = 0. The 
maximum density reached at the beam ener gy o f 50 and 
400 MeV/nucleon is about 1.2po and 2po [33, respec- 
tively. The turnover clearly indicates that the double 
neutron/proton ratio reflects closely the density depen- 
dence of the symmetry energy as shown in Fig. ^ This 
observation also indicates that systematic studies of the 
double neutron/proton ratio over a broad beam energy 
range will be important for mapping out the density de- 
pendence of the symmetry energy. It is useful to men- 
tion that the double 7r~/7r+ ratio in the reactions at 400 
MeV/nucleon has also been examined. There is some 
indication that the double tt^/tt"*" ratio for intermediate 
energy pions is more sensitive to the symmetry energy 
than the single 'k~ j-K^ ratio in the -'^^^Sn-|--'^^''Sn reaction. 
However, unlike the double neutron/proton ratio in the 
same ensemble of events, the double tt^/tt"*" ratio has 
much larger statistical errors mostly because of the poor 
statistics of the single 7r^/7r+ ratio in the ^^^Sn+^^^Sn 
reaction. Studies on the double tt^/tt^ ratio with signif- 
icant more events are in progress. 



IV. SUMMARY 

In summary, within the transport model IBUU04 we 
investigated the double neutron/proton ratio of nucleon 
emissions taken from two reaction systems using three Sn 
isotopes at the beam energy of 50 MeV/nucleon and 400 
MeV/nucleon, respectively. It is found that the double 
neutron/proton ratio has about the same sensitivity to 
the density dependence of symmetry energy as the single 
neutron/proton ratio in the more neutron-rich system of 
the two reactions. Since the double neutron/proton ra- 
tio has the advantage of reducing systematically the in- 
fluence of the Coulomb force and has smaller systematic 
errors, it is more useful than the single neutron/proton 
ratio of nucleon emissions for further constraining the 
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symmetry energy of neutron-rich matter. 
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